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We develop a technique for finding the dynamical evolution in time of an averaged density matrix. 
The result is an equation of evolution that includes an Effective Hamiltonian, as well as decoherence 
terms in Lindblad form. Applying the general equation to harmonic Hamiltonians, we confirm 
a previous formula for the Effective Hamiltonian together with a new decoherence term which 
should in general be included, and whose vanishing provides the criteria for validity of the Effective 
Hamiltonian approach. Finally, we apply the theory to examples of the AC Stark Shift and Three- 
Level Raman Transitions, recovering a new decoherence effect in the latter. 
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I. INTRODUCTION 

The density matrix is the mathematical object that 
carries all measurable statistical information about a 
quantum system, and therefore, completely characterizes 
the state - whether pure or mixed. However, in real- 
ity, the physical perception of any quantum system takes 
place over a finite time interval, rather than instanta- 
neously. Thus understanding the behaviour and evolu- 
tion of the density matrix convolved with some time av- 
eraging function is essential for a complete understand- 
ing of quantum dynamics, as observed in any realistic 
circumstances. 

Moreover, often in practical applications, such as light- 
matter interaction systems, the Hamiltonian has two dis- 
tinct parts - one that oscillates at a high frequency and 
one at a much lower frequency. If one were observing 
the quantum systems with a time resolution which is too 
slow to discern the high frequency effect, one might idly 
suppose that the high frequency component of the Hamil- 
tonian should not play much of a role. And yet effects 
such as the AC Stark Shift or the Lamb Shift, can be 
ascribed to just such high frequency terms in the Hamil- 
tonian. The purpose of this paper is to formalize this 
idea, and find the "Effective Hamiltonian" that has the 
effect of this high frequency component included, and 
examine the validity of this approach in detail. 

The paper is organized as follows. In section II, we 
derive a general formula for the evolution of the time av- 
eraged density matrix, finding the Effective Hamiltonian 
and decoherence terms. Applying this theory to the class 
of Hamiltonians with harmonic disturbances, we derive 
an equation of evolution in the form of Lindblad's open 
system dynamics. The formula for the Effective Hamil- 
tonian found confirms previous results [TJ, and new de- 
coherence terms are found to result from the averaging 
process. 



Finally, we test this theory on known physical systems, 
the AC Stark Shift, and the Three-Level Raman Tran- 
sitions, finding a new decoherence effect in the latter, 
which is potentially realizable through experiment. 

Similar work was completed by the author in the con- 
text of Effective Hamiltonian Theory 1 , and by other 
authors such us Cohen- Tannoudji [2 , and Shore [3J- Ef- 
fective Hamiltonians were also used for various systems, 
such as ion traps and cavities, by James [4], Plenio et al. 
[5], Maimer and S^renson [5J. 



II. EVOLUTION EQUATION OF AVERAGE 
DENSITY MATRIX 



A. The Unitary Evolution 

We start by defining the time-average of an operator 
6(t) as 



6(t) 



f(t - t')6(t')dt' : 



(2.1) 



where f(t) is real- valued, and has unit area. In partic- 
ular, we are interested in the averaged density operator 
p(t). The positivity, unit trace, and hermiticity of p(t) 
imply the same properties for p(t), meaning it too is a 
density matrix. This averaged density operator p(t) can 
be interpreted in two ways: one as a purely theoreti- 
cal construct which we calculate as a mathematical tool 
to discard high frequencies, and the other as physically 
representing the "perceived" density matrix through any 
physical interactions with an apparatus that take a fi- 
nite nonzero time, with the averaging kernel f(t — t') 
representing the strength of this interaction in the time 
window involved. 

By assumption, p{t) represents the evolution of a 
closed system, and therefore is given by the unitary evo- 
lution: 
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p(t) = U{t,to)p(t )U(t,t ) 



(2.2) 



where U(t,to) is the familiar time-ordered evolution op- 
erator which satisfies the Schrodinger equation, viz., 



ih 



dU(t,t ) 

at 



H(t)U(t,t ) 



(2.3) 



In order to calculate a series expansion for U(t,to), we 
adopt the standard approach of replacing H(t) by XH(t), 
where A is a dimensionless expansion parameter. One can 
imagine A gradually being increased from to 1, rep- 
resenting the Hamiltonian being "gradually turned on" . 
When A = 1 we have the standard solution in eq.(2.2|. 
This suggests that we can write U(t,to) as a series in 
powers of A, that is 



U(t,t Q ) = ^2 Xn Un(t,to) 



(2.4) 



Substituting eq.(2.4) into eq.(2.3) and matching the coef- 



ficients of like powers of A, we obtain a recursion relation 
for U n (t,t Q ): 



ih 



dU (t,t ) 

dt 
d U„(t,t ) 

dt 



= F(t)t7„_i(i,i ) 



n > 1. 



(2.5) 



Given that for A — > we must have U(t,to) = I (the 
identity operator), then Uo{t, to) = I. 



Substituting recursively into eq.(2.5l one obtains the 



standard expansion for the time evolution operator. On 
substituting from eq.(2.4) into eq.(2.2), one finds an 
equation of evolution for ~p(t) in terms of U n and powers 
of A. Suppressing for the moment the explicit dependence 
on t and to, and denoting po = p(to), we have 



p = j2\ m+n u mP oui 

m,n 

oo k 



i=o 



J2 xk£ k[p 

k 

£[po], 



(2.6) 



where £fc[po] are linear maps acting on po, and defined as 
the bracketed term in the second line. £[po\ is the linear 
map that gives p. Note that from the definition, £ = X, 
the identity linear map. 



B. Inverse of £ 

To proceed, we would like to invert this transformation 
to find J-\p] = £~ l [p] that operates on p to give po- That 



Po = F[p] 

^A fe J- fc [p]. 



(2.7) 



Naively one might invert the time evolution using the 
unitarity of U(t, to), i.e. by swapping t and to one should 
obtain the inverse. However is this approach is invali- 
dated by the time averaging, since unitarity no longer 
holds. Instead, we use the fact that T and £ together 
give the identity transformation (i.e. J-"[£[p]] = I[p])- 
Thus, postulating that F may be expanded in powers of 
A, and comparing coefficients for powers of the latter, we 
find: 



J2 \ m+n Fm[£n[p]] =A°I[p], 



(2.8) 



m,n— 



which implies 

oo k 

X> fc (£^-;b]])=A :r[p]. (2.9) 

k=o j=a 

Comparing coefficients of powers A, we find the first few 
terms in the expansion of T as follows: 

•Po = £q = 1- 
Ti = -£1 

F 2 = -£ 2 +£ 1 [£ 1 ]. (2.10) 

The recursion relation for higher order terms is provided 
in the appendix. 



C. Approximation up to Second Order 

Next, we find the effective evolution equation of p in 
terms of A, so we can examine the lowest order terms and 
see what they tell us about the Effective Hamiltonian. 
Differentiating eq.(2.6) with respect to time (denoted by 
a dot), and substituting from eq.(2.7) we have 



ih^4^- = ih£[p{t )} = ih£[T[p(t)}] 



Of 



OO K 

= J2 xk {Y, ih ^\.^-Mt)}}}- (2- 11 ) 



3=0 



This can be written as 

.Mt) 



= ^A fc £ fc [p(t)] 



(2.12) 



Evaluating Ek and J~k explicitly, we find, up to second 
order, 



C [p] = iMo^olp}] = 0- 



(2.13) 



where T>2[p] are the decoherence terms, i.e. all the terms 
where p is "sandwiched" in C2' 



V 2 [p] = HpU\ - HpU\ + U lP H + U lP H. (2.22) 



Ci[p] = iHSx^olp}] + iH£o[Fi\p]] 
= Hp- pH, 



(2.14) 



C 2 [p] =HUip-HU lP + Hptj\ - HpU\ 

-pU\H + pUlH-UipH + U lP H. (2.15) 



Eq. (2.12) together with eqs. (2.13) - (2.151 form the 



principal result of this section. Higher orders are included 
in the appendix. The Z\ term is just the Heisenberg 
equation using the averaged Hamiltonian. Note that the 
total expression for Ck is anti-Hcrmitian, as we would ex- 
pect, since the LHS of eq.(2.12 ) is a pure imaginary multi- 



plied by a Hermitian operator, making it anti-Hermitian. 
Also, the familiar pattern of "average of the product mi- 
nus product of the averages" , in a manner reminiscent of 
the definition of covariances, occurs repeatedly. 



D. Comparison with Previous Results 

Previous results due ref. [T], show that to first order, 
the Effective Hamiltonian is given by 



where 



# e// = # + -(! + it), 



A = HU 1 -HU 1 , 



i f = -U 1 H + U 1 H. 



(2.16) 

(2.17) 
(2.18) 



This holds since U{ = —U\. Using the von Neumann 
equation, this implies that to first order: 

■ b d P rTv 1 A + Ai i + lt 

^ = ^ + — ^— P-P— ^— ■ (2-19) 

However, note that the derivation in ref. [I] simply dis- 
carded anti-Hermitian terms to satisfy the Hermiticity 
requirement for the Effective Hamiltonian. Our deriva- 
tion however treats this problem with more rigor, and 
arrives at the detailed solution. From the previous sec- 
tion up to second order (i.e. including C\ and £2), we 
have 

ih^ = [H, p]+Ap- pit + p 2 m f (2.20) 



87i A — /it 

th^ = [H eff ,p] +{——,p}+V 2 \p], (2.21) 



Our results correspond to those of ref. 1 if i = it 
and the decoherence terms are identically zero. As we 
shall see later, it turns out the decoherence terms are in 
general important, and yield an evolution closer to the 
exact case. Thus eqs. (2.12) - (2.15) are the correct ex- 



pressions which should be used in general. The derivation 
in this section relied on the Schrodinger picture, however 
it is valid in the Interaction picture as well. In the latter 
picture however, the interpretation of the density matrix 
changes somewhat. 



III. HARMONIC TIME DEPENDENT 
HAMILTONIAN 

Following ref. pQ , we will now apply our general result 
to a class of harmonic Hamiltonians of the form: 

N 

H{t) = Ho + 2^ h n exp(—iuj n t) + h^exp^ujnt) , (3.23) 

n=l 

where Hq is independent of time. This is representative 
of a wide class of problems, particularly in the interaction 
picture Hamiltonian. In this case, the first term in the 
expansion for U\ is 

Ui(t) - ^j^Ho + Vi(i) - Vr(to), (3.24) 

in 

where 

N 

Vi(t) = ^2 i — lh n exp(-iu> n t) - /^exp(iu; n £) j (3.25) 

n— 1 n 

We assume the averaging kernel /(■) is an ideal low pass 
filter, and an even function in time. We further assume 
that the frequencies ui n are sufficiently high that they 
are filtered out, but sufficiently close to each other so 
that terms oscillating at difference frequencies (uj n — us m ) 
will pass the filter. Thus, we can make the following 
assumptions: 



exp±iw„t = 



exp ±i(uj n + uj m )t = 



exp ±i(ui n — uj m )t — exp ±i(ui n — uj m )t. (3.26) 

This implies that 



H(t) = H Q , 



V 1 (t)=0, 



Ui(t)= { ^-^H Q -mo) 
in 



(3.27) 



Using the Harmonic Hamiltonian eq.(3.23) along with 
eq.(3.26) and our initial results in eq.(2.21 ) and eq.(2.22), 



we obtain the expression 

dp 



ih 



at 



: [h 

N 

El 



tf h 

• ">m""n 



h h) h) h 

P-P 



+ 



huj n huj 

n,m 

HnPhn + hnptlln h^pK 



hui m 
hnphL 



h /i* 
hui n ' 



huj 

H„ff,p 



[ 

oo 1 

N li - 



hu n 



-u„)( 



■•K 



■>} 



•)* n 



h\nP h n + Kph^ l(u 



where 



and 



H, 



eff 



Ho 



1 


1 


i 1 1 
— ± — 

\LU n LO, n 


LU nm 


2 


N 


1 





n<-^nrr 



m i h> n \e 



i{uj n 



i)t 



(3.28) 



(3.29) 



(3.30) 



We can then formulate eq.(3.28) in Lindblad form: 
&p 



th- 



at 



[ H eff,p] 

+ £ T-— ({LlL n ,p} - 2L n pLl 

+ {L n Ll,p}-2LlpL n ), (3.31) 



where 



XJrn, 1 1 j- 



(3.32) 



The Effective Hamiltonian obtained here is exactly the 
same one in ref. pQ. In addition, we have found the ex- 
pression for the decoherence terms. If we have only one 
frequency present in the Hamiltonian, the decoherence 
terms vanish since l/u>~ m = 0, and the Effective Hamil- 
tonian alone perfectly describes the evolution to second 
order. In more general cases, the decoherence term must 
be considered. 

The appearance of decoherence Lindblad terms above 
means the average system can be thought of as an open 
system, even though the underlying "real system" is 
closed, and evolves unitarily. The interpretation of these 
terms in general may be surprising. When one averages 
quantities, one is of course throwing away information, 
and thus increasing the entropy of the system. For ex- 
ample, a low-pass frequency filtering removes information 



about high-frequency processes. Naturally, this increas- 
ing entropy will lead to decoherence terms in the evolu- 
tion equations for the average. 

One point we have neglected so far is the truncation of 
our series at second order. We can obtain a sufficiency 
condition for the legitimacy of this approximation simply 
by considering the ratio of higher order terms of the ex- 
pansion. We accomplish this by assuming that -^- « 1 
Vn, where r\ is the largest eigenvalue of H . This ensures 
higher order terms are progressively smaller, and thus 
may be safely discarded. 



IV. EXAMPLES 

In this section, we revisit two examples cited in ref. 
PQ, namely the AC Stark Shift and Three- level Raman 
Transitions. Note that for both examples, we work in 
the interaction picture [5], starting with the interaction 
Hamiltonian. Inclusion of the decoherence terms given 



in eq.(3.31 1 will illustrate their importance and elucidate 



the validity of the Effective Hamiltonian model. 



A. AC Stark Shift 

Consider a two level atom interacting with an external 
harmonic force. The Interaction Hamiltonian is given by, 



jW*) = ^{|2><l|exp(-»At) 



|l)<2|exp(iAt)}. 

(4.33) 



Applying eq.(3.28l to HAc{t), an d evaluating -ff e //, we 
find 



H, 



eff 



nn 2 

4A 



(|2)(2|-|1)(1|). 



(4.34) 



We also find that all decoherence terms vanish, since 
there is only one harmonic operator, hi = 4p|2)(l|, and 
therefore 1/wj^ vanishes by definition. So in this case, 
the Effective Hamiltonian is time independent. 

To show the Effective Hamiltonian Theory in action, 
we compare both the exact and time-averaged evolution 
of the density matrix for a specific numerical example. 
Recall that we are in the interaction picture, and there- 
fore the density matrix differs from the standard one in 
the Schrodinger picture. The exact evolution is given by 



ih^- = [H AC (t),p(t)}, 



(4.35) 



and time averaged evolution via the Effective Hamilto- 
nian follows 



ih d ^ = [H eff ,p(t)]. 



(4.36) 



We set the ratio b = £ , which represents the strength 
of the applied AC field, and measure time in units of the 



characteristic time ^. Then we start with an arbitrary 
density matrix, and plot the real part of the off diagonals 
- i.e. the real part of the coherences between the popu- 
lations in the two levels. We vary b, and plot the result 
below for b = 0.3. 



Time Evolution of Coherence, b 
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FIG. 1: The evolution of the real part of the off diagonal entry 
for the density matrix, with the interaction strength set at b 
= 0.3. Both axes dimensionless. 



From the close correspondence graphs, it is clear this 
method is an effective first order approximation to the 
frequency and amplitude of oscillation. As expected, the 
time averaged evolution viz the Effective Hamiltonian re- 
sembles the exact evolution without the superimposed 
high frequency components. Note that the relative phase 
of the exact and averaged evolution just depends on ini- 
tial conditions, and zero of time. Initial conditions were 
artificially modified in the above figure so the two graphs 
are in phase. 



B. Three-Level Raman Transitions 

Applying the theory to Raman Transitions, we start 
with the following Interaction Hamiltonian 

H(t) = ^i |3) (I|e- ,;wit + ~ |3) (2|e- ,; " 2t + h.a. (4.37) 




FIG. 2: Illustration of the Three-Level Raman Transitions 



same Effective Hamiltonian 



H eff - 



no 2 

L;i (|3)(3|-|l)(l|)-^(|3)(3|-|2)(2|) 



4w+ 



(|1)(2| 



£(CJl —U)%)t 



|2)(1| 



5 -i(wi-w 2 )q 



(4.38) 



The evolution of the density matrix then follows: 

-&P yfr -1 hVtiVl-2 



ih- 



H, 



ff'Pl 



({|2)(l|,p}-2p 12 |3)(3| 



dt L " J J " J 4w 12 

- 2p 33 |2>(l|y(— »)* - ({|I)(2|,p} - 2p 21 |3)(3| 



2p 33 |l><2|)e 



-i((JJl— U>2)t 



where 



MM 



(4.39) 



(4.40) 



We already assumed that ui\ and (x>i are close in value, 
which means -4r- will be small compared to -ff e // , making 

the decoherence terms relatively small as expected. Now, 
we write p as a generalized Bloch vector, 



p = I + r x X + r y Y + r z Z + r w W 
+ r xa X a + r ya Y a + r xb X b + r yb Y b 



(4.41) 



where 

A = |I)<2| + |2)(1| 

Z=|l)(l|-|2)(2| 

A a = |I)(3| + |3)(l| 
A 6 = |2)(3| + |3)(2| 



y = -*(|i)<2|-|2)<i|) 

TT=^=(|I)(1| + |2)(2|-2|3)(3|) 

y Q = ^(|l)(3|-|2)(3|) 

n = -i(|2)(3|-|3)<2|). 

(4.42) 



The operators above are the Gell-Mann matrices, and 
form a basis for SU(3). The operators X, Y, and Z are 
the standard Pauli operators for the {|1),|2)} subsys- 
tem. Their analogues for the {|1), |3)} and {|2), |3)} sub- 
systems are labeled by the subscripts a and b respec- 



tively, and will soon be discarded. Substituting eq.(4.41 ) 
in eq.(4.39), it simplifies to the following equation: 



dp ilin 2 ' 



Z(r x sin 9 + r y sin 9) — r z (Y cos 9 + X sin 9) 



-(^-^)(r x Y-r y X) 
4 v wi uj 2 



V3 



j^i0 2 r 



2w 12 
+ W(r x sin 



,(Xsm9 + Ycos9) 



(4.43) 



where 



Applying our formula eq.(3.28l, we once again get the 



(wi -uz)t, 



(4.44) 



and g is some linear functional that only depends on 
r xa ,X a ,r ya ,Y a ,r xb ,X b ,r yb , and Y b . This implies that 
the coherences between level 3 and and the other two 
levels (i.e. the (1,3), (2,3), (3,1), and (3,2) entries of 
the density operator) form a closed subsystem, and only 
affect each other, evolving independently of the rest of 
the density matrix. Keeping this in mind, we ignore this 
subsystem entirely, and focus on the evolution of r x , r v , 
fz, r wi which together form a 4 dimensional analogue of 
the Bloch vector. This vector follows the evolution of the 
following dynamical system, viz: 



Or 

at 



= Av 



(4.45) 



where 



.4 






—a 


-/3sin0 


—7 sin 


Q 





—j3 cos 9 


—7 cos 


3sin# 


(3 cos 9 








7sin# 


— jcosO 









(4.46) 



The constants, a, /3, 7 are given by 



$■■■ 



!'■ 



4 y uji 

n x n 2 

2lo+ 



Of 

L0% 



V3 



2 
2u7n 



(4.47) 
(4.48) 
(4.49) 



As we shall see, we can view a and /3 effectively as driving 
terms, where 7 represents reduction in frequency due to 
decoherence. To simplify the matrix A above, we go to 
the co-rotating frame, and define: 



where 



M s = 



f = MgY, 



cos 9 -sin 9 

sin6> cos6» 

10 

1 



(4.50) 



(4.51) 



We can then define the vector d as the first three com- 
ponents of r, leaving f w as a scalar. We also define new 
vectors 7 and $7 as follows 



7 



n = 



(3 


■Wl 



Our equations then reduce to the simply form 

d = n x d — r w ~f 

f w = — V d. 



L0 2 

(4.52) 



(4.53) 



This presentation intuitively describes the evolution of 
the system as the rotation of the bloch vector about Q 
through the term given by the cross product (as is al- 
ready well known with the Rabi oscillation [TU]) plus a 
perturbation given by the terms in 7. Now to solve this, 
we define the unit vectors en and e 7 as the orthogo- 
nal unit vectors along $1 and 7 respectively. Together 
with e p = eoxe 7 , they gave an orthonormal basis for 3d 
space. Using this basis, we write: 



d 

n 

7 



dnen 

7e 7 



ti-y tJ/y 



dpQp 



(4.54) 



substituting this ansatz for d in eq.(4.53) and resolving 



the equation along each unit vector, we get the following 
system of coupled differential equations: 






ir u 



-7dv, 



from which we get 



d~, 



-{VL 2 



7 2 K 



(4.55) 



(4.56) 



Solving this simple second order ordinary differential 
equation, we have 



cL = R cos tot 

n 

d p = R— sm cat 

LJ 



f w {i) 



7 
-it— sin tot 



7 - 

fwO 



(4.57) 



where R is some constant oscillation amplitude, and 
r w0 = r w (0) is a constant, both dependent on the ini- 
tial conditions. The frequency of oscillation lo given by 



= n 2 - 7 2 

= (a + LOi -w 2 ) 2 +/3 2 -7 2 



(4.58) 



Note that we set the zero of time such that there is no 
phase angle in the argument of the trigonometric terms. 
So our final solution is 



d(i) = en<fc 



„ 7~ 



i?(e 7 cos cot 



e„— sm Lot) 

CO 



(4.59) 



Firstly note that when 7 vanishes, we have lo = Cl, and 
Bloch vector oscillates at the Rabi frequency. When 
f-wo = 0, this oscillation is around the torque vector 17, 
as in the well known result [ID] . As 7 increases, it re- 
duces the frequency of precession of the Bloch vector, 
and changes the oscillation path from a circular one to 
an Elliptical one. As f w o increases, the centre of this 
oscillation shifts in the e p a direction (perpendicular to 



both $7 and 7. The elliptical nature of the oscillation 
means the length of the vector is not constant, rather it 
oscillates at frequency u. Note that the length squared 
of the Bloch vector corresponds to the trace of the square 
of the underlying density matrix - i.e. a measure of its 
purity/coherence. If we define I 2 as 



l 2 = d 2 a + £ f + d 2 p . 



(4.60) 



Taking its time derivative and using [4~55} we find 



d(l 2 ) 

dt 



2(rfndo - 
= -2jf w d 



LL-y (X^-y 



dp&p) 



7 2 2 

— R sin 2ujt — 2 r )Rf W Q cos ust 



(4.61) 



So in general, the length of the Bloch vector (i.e. the 
coherence of the effective 2 level system) oscillates at fre- 
quency u, and if r w o vanishes, it oscillates at 2w. 



V. CONCLUSION 

We have demonstrated a more rigorous derivation of 
Effective Hamiltonian Theory and its domain of applica- 
bility. Additionally, it has been shown that additional 
terms are created resembling Lindblad evolution for har- 
monic Hamiltonians - implying that the averaging pro- 
cess introduced a small decoherence factor. Applying 
this theory to examples such as the AC Stark Shift and 
Raman Transitions, we find it introduces some minor cor- 
rections. 

In the future, applying this theory to other known sys- 
tems and testing its limitations would inform us about its 
usefulness as a tool. In particular, applying it to systems 
where the newfound decoherence terms play a larger role 
would help us better understand their exact role. Ad- 
ditionally, it would be useful to come up with a general 
interpretation of the Lindblad terms. For example, the 
averaging process over a given system can be seen as ob- 
serving an analogous system-reservior pair for some hy- 
pothetical reservoir, and then observing the system alone 
while "tracing out" the reservoir. 



Acknowledgement 

This work was funded by the Natural Sciences and 
Engineering Research Council of Canada (NSERC), and 
the University of Toronto Physics Department. 



Appendix A: Derivation of Tu and Lu terms 

The recursion relation for calculating J- n is given by: 



ra-l 

Fn = ~Y, J: A E n-Ap]]- (1-1) 

3=0 

In the following, will calculate £& and T^ terms for the 
first few orders. Using eq.(2.6l and we get 



£o[p] = P 
8 1 [p]=W 1 p + pWl 

£2[p] = &2P + UipUl+pU 1 2 



£a\p] = U aP + U 2 pU\ + tJ x pU\ + pU, 

Applying the differential operator ihjL: 

iHSq [p] = 

ih£\ [p] = Hp — pH 



(1.2) 



ih£ 2 [p] = HUip + HpU{ -U x pH- P U\H 



ih£ 3 [p] = HU 2 p + HU lP U\ - U 2 pH + Hp(j\ 



- UiplllH - pU\il. 



Finally from eq.(2.7l 



Mp] 



(1.3) 



(1.4) 



T 1 [p] = -U 1 p-pUl 

t 2 [ p ] = -¥ 2 p - Uxpul - piil + ih(ihp + p&r 

%p + pUl)ul. (1.5) 



The third order term £3 is found to be 

£3 [p] =HU 2 p~HU 2 p + HpUl - HpUl 
- pU\H + pUl H - U 2 pH + U 2 pH 



MipU\ - HUipUl - HUipU\ 



HU-LPUI + 2HU-LPUI 



UipU\H + UxpU\H + U lP U\H 
UxpUlH - 2U lP UlH 



HU X Uxp + HUx UxP - HpU\ul + HpU\ Jj\ 



+ pU\ U\H - pU\ U{ H + UxlIxpH - Ux UxpH. 

(1.6) 
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